Phase transition of the nucleon-antinucleon plasma at different ratios 
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We investigate phase transitions for the Walecka model at very high temperatures. As is well 
known, depending on the parametrization of this model and for the particular case of a zero chemical 
potential ( ^ ), a first order phase transition is possible [Tj. We investigate this model for the case in 
which n ^ . It turns out that, in this situation, phases with different values of antinucleon-nucleon 
ratios and net baryon densities may coexist. We present the temperature versus antinucleon-nucleon 
ratio as well as the temperature versus the net baryon density for the coexistence region. The 
temperature versus chemical potential phase diagram is also presented. 



00 

o 

O 
(N 

Oh 
D 



I 



> 

in 

(N 

On 
O 
00 
O 

> 

X 

u 



PACS number(s): 21.65.+f,24.10.Jv,21.30.-x,21.60.-n 

Finite nuclei and nuclear matter properties have been 
reasonably well described by the Walecka model |2j . This 
renormalizable model employs nucleons and mesons (ct 
and w) as the degrees of freedom. The sources for the 
fields are the scalar (p^) and vector (p) densities associ- 
ated with the Lorentz scalar (S) and vector {V) interac- 
tions. It is well known that this model, after the fitting 
of the experimental infinite nuclear matter binding en- 
ergy at the saturation density, gives for this system a 
high bulk incompressibility, as well as a high spin-orbit 
splitting energy, when applied for finite nuclei. In order 
to overcome this problem, nonlinear relativistic hadronic 
models which include cubic and quartic scalar mesonic 
sef-couplings have been proposed [3]. The preliminary 
success of this kind of nonlinear Walecka model, estimu- 
lated the proposal of different parametrizations to bet- 
ter describe nuclear matter properties [4|. Still aim- 
ing at refining its predictions for nuclear matter and 
neutron stars properties, new kinds of models including 
scalar (vector) self-coupling interactions added to usual 
nonlinear Walecka models [31 [4] have been constructed 
[5]. Effective hadronic models with density dependent 
coupling constants have also been proposed [3]. 

In general, (non)linear Walecka models models 
m [6] have been investigated under extreme density and 
temperature regimes. Most of such studies focus on the 
hadronic quark-gluon plasma phase transition, which is 
not the aim of our present work. Here, our purpose is 
to extend the very interesting work of Theis et al [1] . In 
their study T , the Walecka model was investigated in the 
extreme high temperature regime in which the chemical 
potential ( )is zero. This leads to a situation in which 
the net baryon density ( p ) also vanishes. In this case, the 
number of nucleon and anti-nucleon are the same. This 
study showed that, in this nucleon-antinucleon plasma at 
very high temperature, a phase transition exists. 

In this work, we will investigate the Walecka model 
under an extreme temperature regime, but allowing the 
number of anti-particles and particles to be different in 
the two phases it can lead to. The question we pose here 
is whether by allowing different ratios of nucleons-anti- 
nucleons a phase transition scenario is still present, as in 



the particular case p = 0, still in the so called 'no-sea' 
approximation. That is, we consider explicitly only the 
valence Fermi and Dirac sea states. Let us here remark 
that such an approximation, as investigated carefully by 
the Ohio group [7 , becomes good at low energies, where 
the contributions of the Dirac sea can be renormalized in 
effective coupling constants. Of course, a rigorous way 
to perform our investigation would be to include vacuum 
polarization effects in the Walecka model for high exci- 
tation energies, where the anti-particles play an essential 
role. This treatment, however, is beyond the scope of 
this work and we perform an exploratory investigation 
even though the validity of the 'no-sea' approximation in 
this regime may not be granted. 

Nowadays, high energy experiments reveal evidence of 
nuclear systems with very small baryonic density in the 
study of particle yields measured in central Au — Au col- 
lision at RHIC. Different experimental analysis (STAR, 
PHENIX, PHOBOS, BRAHMS) furnish the antiproton- 
proton ratio p/p w 0.65 for a temperature of YlAMeV 
and a chemical baryonic potential of AQMeV estimated 
from thermal models to fit antiparticle-particle ratios [8] . 
High energy Ph — Pb collisions show that this ratio reach 
values around 0.9 If relativistic hadronic models are 
to be used in the description of the multiplicity observed 
today in ultra-relativistic heavy ion collisions [TOj ITT| . 
the behavior of these models at high temperature regimes 
may be of importance. 

The thermodynamics of the Walecka model may be 
given, for finite temperature, in terms of the energy den- 
sity and pressure functionals. 
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and In the expressions above, 7 is the degeneracy factor 
(7 = 4 for nuclear matter and 7 = 2 for neutron mat- 
ter ), rife and stand for the Fermi-Dirac distribution 
for baryons and antibaryons respectively, with arguments 
{E*^v)/T . E*{k) is given by E*{k) ^ (P + M*2)i/2 ^ 
whereas an effective chemical potential, which preserves 
the number of baryons and antibaryons in the ensemble, 
is defined by u = 11 — V , V ^ C'ljPb/M'^, where ^ is the 
thermodynamic chemical potential. The solution for the 
equation of state is obtained explicitly through the mini- 
mization of £ relative to the scalar field, or equivalently 
to m* = M*/M, 
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This equation, known as the gap equation, has to be 
solved self-consistently with Eqs. [l]|2] and provides the 
basis for obtaining all thermodynamic quantities in the 
mean field approach we are using. 

Usually, the constants and are determined in 
favor of the experimental nuclear matter binding en- 
ergy {16MeV) at the saturation density {ph = Po = 
0.15fm~^), at T = 0. At finite temperature, this model 
predicts a critical liquid-gas behavior as a van der Waals 
EoS. Its critical temperature is around 18MeV [2]. For 
this temperature, the antinucleon-nucleon ratio, given by 
R = p/ p is negligible. Only when the temperature gets 
higher, the antinucleon density starts to take significant 
values. 

In a very interesting work [T] , the Walecka model was 
studied in the extreme situation in which R — 1 . In this 
case, Pb = and v ^ p. ~ 0. This study showed that, 
in this nucleon-antinucleon plasma at very high temper- 
ature, a phase transition exists. It is also remarkable 
to note that the order of the phase transition itself be- 
comes dependent on the versus space parameter. 
By small changes on these parameters (which means to 
change by a few percent the nuclear matter binding en- 
ergy and the saturation density), the phase transition 
changes from first to second order. If we take the val- 
ues of Cl = 359.35 and = 275.12, fitting the infinite 
nuclear matter binding energy as 16MeV at a density 
of 0.15/m~'^, the phase transition is of first order. The 
findings of Ref. [T] may be summarized in Fig. [l] 

In Fig. [T| we see that the nucleon effective mass de- 
creases abruptly at T « ISiMeV. Since all other ther- 
modynamic quantities are dependent on to*, this effect 
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FIG. 1: The effective mass, entropy, energy density, pressure, 
and the specific heat for i? = 1 as a function of the tem- 
perature. The mass is in units of the nucleon mass and the 
entropy, energy, and pressure are in Stefan-Boltzmann units. 



also manifests itself in energy density, pressure, and spe- 
cific heat, calculated as the temperature first derivative 
of the energy density. The fast increase of the entropy 
is a clear signal of a first order phase transition. Since 
the baryonic density is kept constant pi, — 0, the order 
parameter should be the entropy. It is a very curious 
constrained system, interpreted as a nucleon-antinucleon 
plasma |T|. 

Here, following pj, we start to study the Walecka 
model in different constrained cases. The main question 
is to understand the behavior of some thermodynamical 
quantities when the ratio R = p/p varies. Now, all 
quantities will have contribution from the baryonic den- 
sity, contrary to Fig. [T| where pb = 0. 

Let us now remark that we have fixed numerically 
the ratio within a precision of one part in a thousand. 
It means that the set of equations (1-7) is solved self- 
consistently with R ~ i?/1000 < R< R + i?/1000. In 
figure [2] we present M* /M as a function of the temper- 
ature for different values of R. 

For the same values of R investigated in Fig. |2] the en- 
tropy behavior as a function of the temperature is shown 
in Fig. |3] 

Both the nucleon effective mass and the entropy fol- 
low the same abrupt decreasing (increasing) behavior one 
sees for the particular case R — 1. From these figures 
we see that, as the ratio p/p decreases, the temperature 
in which an abrupt behavior arises also decreases, but 
keeping the same character. In principle, there is an im- 
portant difference between the system described by Fig. 
[1] and those of Figs. [2] and |3] In the first, R = \ and 
Pb — Q along the temperature variation. In the second, 
the ratios are kept constant while pb varies. Visually, 
Figs. 2-3 suggest phase coexistence at some tempera- 
ture range. According to the Gibbs criteria, if one has 
phases 1 and 2, a phase coexistence arises when pi = P2, 
pi = p2 and Ti = T2. The critical temperature is 
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FIG. 2: The effective mass as a function of the temperature 
for several values of R. From the left to the right: R = 
0.1,0.3,0.5,0.7,0.9. 
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FIG. 3: The entropy, in Stefan- Boltzmann units, as a function 
of the temperature for several values of R. The values of R 
are the same of figure [2] 

achieved when, above that, no more phase coexistence is 
possible. For the case = I, a phase coexistence ex- 
ists at T = 183.25MeF and = . It means that, if 
the temperature increases or decreases from this value, 
the phase coexistence disappears as we can see in Fig. 
|4] where the minima of the thermodynamical potential 
have different values for the same temperature. 

This kind of investigation was performed, for example, 
by Asakawa and Yazaki I13| when working with the NJL 
model at finite temperature, where a phase transition is 
also present. Following their procedure, we can verify 
whether the system exhibits phase coexistence or not. 

Surprisingly, despite the signals of phase coexistence 
shown by Figs. [2][3] we did not find phase coexistence 
for a single value of i? when i? ^ 1 . In this case, when 
the system was in thermal (Ti = T^) and mechanical 
(Pi = P2) equilibrium, chemical equilibrium (^1 = /i2) 
was not found. On the other hand, along the apparent 



FIG. 4: Thermodynamical potential, in units of the Stefan- 
Boltzmann pressure, as a function of the nucleon effective 
mass, for temperatures around the critical temperature in the 
case where = {R = 1). 

phase coexistence region signalized in Figs. |2][3]by three 
different values of M* /M and entropy for the same tem- 
perature, we found a stable phase (global minimum of the 
thermodynamical potential) . It indicates the presence of 
stable phases up to the start of the backbending of the 
curves, and new stable phases by any increase of the tem- 
perature, but with a dramatic decrease (increase) in the 
effective mass (entropy). Such a phase transition with- 
out phase coexistence is rare in physics. Therefore, we 
decided to leave the fixed ratio scenario to analyse what 
happens if the system becomes free of such constraint. 

Next, we proceed to investigate the phase coexistence 
in the Walecka model still at high temperature but with- 
out any fixed R constraint. In Fig. [5] we display the ther- 
modynamical potential density —p and f{M* /M) given 
by Eq. ^ and Eq. ([7| respectively. 

In the curves presented, the Gibbs criteria Ti = T2, 
Pi = P2) and /zi = /i2 are clearly satisfied. As we 
continue to decrease the temperature, the Gibbs crite- 
ria can be fulfilled since the chemical potential fi in- 
creases. The curves for the thermodynamical poten- 
tial for T < 182.9Mey (not shown) becomes flat- 
ter as the temperature decreases. This happens until 
T w ISOMeV and ^ w 274:MeV . Therefore, the coex- 
istence region for the Walecka model, with the coupling 
constants given previously and without any R constraint, 
is 180.1 < T < 183.25Mey with < ^ < 274:MeV . 
The phase diagram T x /i is given in Fig. |6] 

Now, we have the values of M* /M which allow the 
phase coexistence for the phase diagram of Fig. [6] With 
these, we can extract the net baryon densities p as a 
function of T which allows phase coexistence. The re- 
sults are presented in Fig. [7] 

Following the same procedure, in Fig. |8]we show the 
different ratios R for which the system affords coexis- 
tence. It is interesting to observe that no coexistence 
exists if one of the ratio is not greater than 1/2. As we 
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FIG. 5: The same as in the previous figure, but showing the 
phase coexistence for temperatures between T = 183.25AfeV 
(upper curve in top panel) and T = 182.9071^6^ (lower curve 
in botton panel). The values of /i are those that make the two 
minima to be have the same value for each temperature. Note 
that, now, R is not constrained and have a different value for 
each temperature and chemical potential. 



FIG. 7: Temperature versus net baryon density for the phase 
coexistence region. 
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FIG. 6: Temperature versus baryon chemical potential for the 
coexistence phase region. 



can see, only in the particular case i? = 1 , which corre- 
sponds to ph = /i, there is a phase coexistence allowing 
only one fixed ratio. 

In short, we start by studying the Walecka model at 
high temperature regime constraining the anti-nucleon- 
nucleon ratio (i?) to be constant but not equal to one as 
has been done by Theis et al. |1, . We have seen that the 
visual signals of phase transition for the effective nucleon 
mass (M*) and entropy (5) versus temperature (T) for 
the cases 1 are tipically the same of those obtained 
in the case R — 1. By this we mean an abrupt decrease 
(increase) of both (M* and S ) iov T > ISOMeF. Sur- 
prisingly, however, we could not find for R ^ 1, contrary 
to the R = I case, a phase coexistence signature by using 
the Gibbs criteria. Phase transition without phase coex- 
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FIG. 8: Antinucleon-nucleon ratio (R) versus temperature for 
the phase coexistence region. 



istence is rare to occur in physics. As far as we know, 
they are theoretically predicted for the anomalous two- 
dimension Kosterlitz-Thouless phase transition [TT, as 
an example. Since we could not find a clear signature for 
phase coexistence at fixed i? ^ 1 we have investigated 
phase coexistence for non-constrained values of R. Our 
study shows that the only phase coexistence regime for 
a fixed R occurs at i? = 1 or, what is the same, /i = 0. 
In the the Walecka model parametrization we have used, 
coexistence phases occur approximately in the interval 
180 < T < 183.25Mey and different values of R at each 
phase are needed. 

If one uses hadronic models (TUl E] to study high en- 
ergy heavy-ion collisions, the models face low density 
and high temperature regime. In the chemical freeze- 
out, anti-nucleon- nucleon are produced at ratios which 
depends on the center-of-mass energy [12] . Although the 
Walecka model is to simple to deal with, it may anticipate 
roughly what can happen with more realistic hadronic 
models. Theis et al [T] analysed (/j, = 0) that, depending 
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on its parametrization, the Walecka model shows first or 
second order phase transition. The same happens with 
different hadronic models [TS]. Therefore, it is to be ex- 
pected that models (with /i = 0) which present first order 
phase transition would follow a behavior similar to the 
Walecka model (Figs. [6]-[8| regarding phase coexistence 
regions. Therefore, depending on the values of T and /i 
used to fit the freeze-out data, the chosen hadronic model 
may be dealing with hadronic phase coexistence region. 



as depicted in figures 6]|8 
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